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ABSTRACT 


Insulation system in high voltage power equipment plays an important role for 
the reliability of the electric power system. So, it is necessary to assess its 
performance to prevent any sudden interruption in the power system. Partial 
discharge (PD) in solid dielectrics could occur due to the presence of a cavity 
or a crack within the insulating material which could be formed during 
manufacturing, installing or/and operating conditions. Since internal cavities 
are the main source of the PD activities, they can lead to causing deterioration 
of the insulation system and consequently a complete failure may be occurred. 
In this paper, PD measurement and simulation are performed on a rubber 
insulating material with variable cylindrical cavity diameters. High-frequency 
current transformer (HFCT) technique is introduced for PD measurement. 
Therefore, a PD simulation model is developed using COMSOL Multiphysics 
program interlinked with Matlab software in order to investigate the influence 
of cavity geometry on the PD behavior in insulating material. Both measured 


and simulated results indicate that PD magnitude is strongly depending on 
changing the cavity size inside the insulating material. The results show a good 
agreement between experimental and simulated data outputs in terms of 
maximum PD magnitude. 
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1. INTRODUCTION 

Partial discharge (PD) is an internal or external electric spark occurrence that does not pass the 
distance among electrodes inside an electrical insulating material during electric field strength [1]. PD usually 
happens due to several defects such as internal small voids, cavities, or cracks inside the insulating material 
formed during the manufacturing process or during the operating service under electrical stress for a long term. 
Under operating conditions, the electric field inside the gaseous cavities exceeds the field of the insulating 
material and this returns to the higher permittivity of that dielectric material [2]. PD may occur if the electric 
field of the gaseous air in the void is enough to be large and the strength of the gas inside the void is 
increased [2]-[4]. During PD occurrence, the gas-filled cavity converts its insulation property to a conducted 
material causing decreasing the electric field inside the cavity to a lower value in a very small duration of 
time [5]. Since internal voids or cavities are the main sources of PD events, insulation degradation in the system 
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could be happened and hence, a complete failure may be occurred [6]. As a result, monitoring of PD activities 
is very significant tool for evaluating the performance of the insulating system. 

Several experimental detection methods have been used on the basis of both electrical and 
non-electrical techniques for detection, analyzing, and estimation of PD activities in solid dielectrics as 
described in [7]-[13]. On the other hand, there are PD models that have been introduced to study the PD 
performance inside cavities in insulating materials. The abc-model is a traditional model that was used for 
modeling and simulation of PD mechanism in solid dielectrics. In this model, the insulating material and the 
cavity are represented as capacitors [8]. Although this model has been used to clarify some observed 
experimental data, it is not realistic in describing the cavity discharge process [14], [15]. Recently, with 
emerging modern technology and upgrading of PD software models, finite element analysis (FEA) software 
programs has been developed to simulate the PD behavior inside the dielectric materials and to understand the 
PD phenomena more deeply. This advanced tool is widely used for modeling of PD based on the geometry and 
the size of the cavity in solid dielectrics [1], [2], [6], [14], [16]-[19]. The model introduced in [6] studied the 
PD events within a spherical cavity at variable applied voltages and frequencies and the simulated results were 
compared to the experimental data. In addition, the model developed in [1] studied the PD behavior on a 
cylindrical cavity within polycarbonate plates. However, there is a shortage of usage of this advanced model 
on some applications such as outdoor insulators due to their complicated geometry as some researchers have 
studied only the distribution of potential and electric field on high voltage outdoor insulators [20], [21]. 

In this work, an experimental study of PD detection and measurement is carried out on a rubber solid 
dielectric material with variable artificial cylindrical cavity diameters. Also, COMSOL Multiphysics is used 
as a FEA tool which is interfaced with MATLAB program software to simulate and model the PD events 
resulting from the cavities inside the rubber insulating material. In addition to that, a comparison study is 
introduced between the measured and simulated results. 


2. EXPERIMENTAL SETUP AND RESULTS 
2.1. Test specimen 

The samples under test were prepared with 5 cm diameter and 5 mm thickness. In order to obtain PD 
sources in the test samples, artificial cylindrical cavities were introduced in the rubber samples with variable 
diameters as 2 mm, 4 mm, 6 mm, and 8 mm in addition to the healthy one. To create the test object, three 
samples are used to be compressed together between two electrodes; the first one is connected to the upper 
cylindrical electrode, the second one is connected to the lower cylindrical electrode, and the third one 
containing the cavity is placed between the two other samples as shown in Figure 1. The upper and lower 
electrodes are cast in epoxy to prevent any discharges between electrodes and the samples. So, the test object 
which contains the three samples represents a rubber insulating material with an internal void or cavity. During 
the test, the test object with electrodes are placed in a polymer holder to be compressed together and to avoid 
any separations between the samples. 
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Figure 1. Test object with cylindrical cavity: (a) schematic diagram of PD model and (b) experimental test 
object 








2.2. Test setup 

A high voltage dielectric test set with 100 kV step up transformer, 50 Hz, (phenix technologies) is 
used as AC source for applying voltage on the test object. The higher electrode is interlinked to the high voltage 
terminal of the transformer and the lower one is grounded. PD is detected by High-frequency current 
transformer (HFCT) (CT100) which is placed around the ground terminal and connected to MPD540 
acquisition unit which is connected to a battery. A fibre optical cable is used to connect the MPD540 unit to 
the USB502 then to a personal computer as shown in Figure 2. 

An alternating current (AC) voltage, 50 Hz, is applied to the test object to produce PD activities 
resulting from artificial cavities. HFCT operation depends on the detection of current pulses with high 
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frequency that could be happened during PD in the test object. These current pulses are converted into charges 
which usually expressed in pico coulombs (pC) [22], [23]. 


2.3. Results 

According to the schematic diagram shown in Figure 2, the apparent charge is measured for each 
cavity diameter in addition to the healthy sample at laboratory ambient conditions. Before performing the test, 
PD measuring system shall be calibrated by a reference PD calibrator (LDC-5/R) as recommended in IEC 
standards and corresponding researches [22], [24]. Figure 3 shows the example of PD level measurements 
pattern with artificial cavities for diameter 4- and 8-mm. Apparent charge is measured with different applied 
voltages for each test object as shown in Figure 4. It has been seen that the PD magnitude increases with 
increasing both applied voltage and the cavity diameter. The charge value of the sample without cavity is the 
lowest value and the sample with 8 mm cavity diameter is the highest value. Also, it has been observed from 
this experiment that the inception voltage is noticed to be around 5.45 kV, 5.23 kV, 5.01 kV, and 4.82 kV for 
cavity diameters 2, 4, 6 and 8 mm respectively. 
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Figure 3. Example of PD measurements pattern at 30 kV, for: (a) 4 mm and (b) 8 mm 
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Figure 4. PD magnitude at variable applied voltages for all samples 


3. PD MODEL SIMULATION 

In this study, a PD model is developed using COMSOL Multiphysics software as a FEA tool which 
interfaced with MATLAB software. The geometry of the model is carried out in two dimensional (2D) axial 
symmetric depending on the unit under test that was used in the measurement section with different cavity 
diameters. The model consists of four main parts; a homogenous insulation material (rubber), a cylindrical 
cavity with 2, 4, 6 and 8 mm diameters, a cavity surface of 0.1 mm thickness, and two copper electrodes as 
shown in Figure 5. 

After adding the data for all used materials from COMSOL library or by manual, boundaries for each 
section are created. The upper copper electrode is set as boundary potential (V= V app. Sin (1 00zt)) and the lower 
electrode is set as boundary ground (V = 0). After that, the meshes in the model are refined to obtain more 
accurate results [2]. 

In this model, it is assumed that the insulation material (rubber) is homogeneous and the PD events 
occur at the cavity center due to the strongest concentration of the electric field in the cavity center [2], [16]. 
The electric potential and electric field distributions through the developed model are solved by (1) [1], [16], 
[19] by choosing ‘electric current (ec)’ from ‘AC/DC module’. For sinusoidal applied voltage, time dependent 
study is selected to be used for solver settings. 


y: |-ovv = 2 (€,2,VV)| = 0 (1) 


Where: V is the electric voltage, o is the conductivity, €, is the vacuum permittivity, and €, is the relative 
permittivity. 
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Figure 5. 2D geometry of the model in COMSOL program 
(1: insulation material, 2: cavity surface, 3: cavity, 4: potential electrode, 5: grounded electrode) 


3.1. PD occurrence 

In this model, two necessary conditions are required for PD occurrence in a cavity [1], [2], [14], [16], 
[17], [19], [25]; the first condition is the electric voltage in the void shall be greater than the intensity of the 
breakdown of the gas which corresponds to the cavity inception voltage (Vinc), and the second condition is that 
free electrons must be available enough inside the cavity to start an electron avalanche process. So, the first 
condition can be achieved by testing that the cavity voltage (Vcay) is greater than the voltage inception (Vinc) 
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that required to initiate PD event. For the second condition, free electrons in the cavity in this study are obtained 
only from cavity surface emission as in [1], [16] depending on the electron generation rate (Ne) according to 
(2) [1], [16], [19]. 


Ne(t) = Negexp( ean (t) /Vine (O) 1) (2) 


Where: Veay is the voltage through the cavity center, Vinc is the inception voltage to initiate PD event, and Neo 
is a constant. 

In (2) is very simple for electron generation rate based on Ritchardson-Schottky law [19], [26]. 
Sometimes the number of free electrons is poor for starting electron avalanche. So, the electron generation rate 
and PD occurrence are considered as a stochastic process [1]. Consequently, the probability (P) of generating 
a free electron in the cavity in a time step (dt) is assumed to be N,(t)dt according to (3) [1], [16], [19]. Then the 
acquired probability value is compared with a random number (R) which has a uniform distribution from 0 to 
1 (O <R <1). At each time step there is a possibility of PD occurrence if V cav > Vinc and P>R. 


P = Ne(t) (3) 


Since the time of the discharge is very short almost in a little nanosecond, the electric potential in the 
cavity drops sharply within a short time period [16]. After a PD event, the electric potential inside the cavity is 
decreased causing the released electrons to lose their accelerated energy and the conductivity and the current 
decrease inside the cavity. When the electric voltage inside the cavity reduced to become lower than the voltage 
extinction (Vext), the PD activity stopped [2]. 

In this current study, PD activity inside the void is modeled dynamically with increasing the 
conductivity of the cylindrical cavity according to (4) [1], [19]. 


i | ae during discharge (4) 
Oo no discharge 


Vcav 
Vinc 








I 
+ eax) 
Io 


Where: do is the initial cavity conductivity at no discharge, Ieav 1s the current through the void, and Io is the 
initiation current at an electron avalanche onset 
The current in the cavity during discharge can be calculated according to (5) [1], [19]. 


(5) 


leav 


f Jf J.ds during discharge 
0 no discharge 


Where: J is the current density through the cavity. Finally, the physical charge value (q) due to PD event is 
calculated according to (6) [2], [19]. 


fi le ae (6) 


4. SIMULATION RESULTS AND DISCUSSION 

The parameters that used in this model are assumed as shown and described in Table 1 represented in 
a MATLAB interactive menu at the beginning run of the model. The inception voltage (Vinc), the extinction 
voltage (Vext), and the initial electron generation rate (Neo) are selected to be compatible with the measured 
results at 50 Hz. Figures 6(a) and 6(b) show the distribution of the electric potential and electric field in the 
rubber material and the cavity in case of 8 mm diameter, for example, before the occurrence of PD in the cavity 
respectively. It has been observed from Figure 6(b) that the electric field inside the cavity is larger than that in 
the surrounding rubber material because of the higher dielectric constant of the rubber material. Therefore, the 
electric field on the top and bottom of the surface wall of the cavity nearest to the electrodes is low because the 
applied electric field is perpendicular to the surface of the cavity [2]. 

Furthermore, Figure 7(a) shows the electric potential after PD occurrence. Also, Figure 7(b) shows 
the electric field distribution after PD event. It has been observed from Figure 7(b) that the electric field across 
the cavity seems to be less than that of the surrounding rubber material because of the dynamic motion of 
electrons inside the cavity [2]. In addition, the electric field on the top and bottom of the cavity surface closest 
to the electrodes is high. After PD occurrence, charges accumulate on the surface wall of the cavity producing 
an opposite electric field which decreases the electric field inside the cavity [2]. 
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Table 1. Definition of used parameters in the model 


Symbol Value Description 
Vap 30 kV Test voltage 
f 50 Hz Rated frequency 
V ext 2 kV Extinction voltage 
Eiiat 4 Permittivity of rubber material 
Eese 1 Cavity permittivity 
Es 4 Cavity surface permittivity 
Ecou 1 Copper permittivity 
Sag 1x101? S/m Rubber conductivity 
Saw 1x10! S/m Cavity conductivity 
Ss 1x10" S/m Cavity surface conductivity 
Sa 5.813x10’ S/m Copper conductivity 
Scie 1x107 S/m Max. conductivity of cavity 
Neo 650 I/s Initial electron generation rate 
Ip 1x10? A Initiating current at an electron avalanche onset 
Contour: Electric potential (kV) Surface: Electric Field Norm (kV/mm) 
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Figure 6. Distributions before the PD event, for (a) Electric potential and (b) Electric field 
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Figure 7. Distributions after the PD event, for (a) Electric potential and (b) Electric field 


Figures 8-11 show the simulation of Vapp, Vcav, Vinc, Vext, and charge magnitude against the phase 
angle during one cycle for 2 mm, 4 mm, 6 mm, and 8 mm cavity diameters. Initially, the cavity voltage starts 
from zero with the applied voltage, and then it exceeds as the applied voltage increases. When the cavity voltage 
arrives to the inception voltage and more released electrons are available radiated from the surface of the cavity, 
the PD can occur. Immediately after PD occurrence, the cavity voltage drops sharply until reaching to the 
extinction voltage where the discharge stops. Increasing the cavity conductivity causes the voltage inside the 
cavity again to be increased and the discharges repeated. It has been seen from these figures that, PD occurs 
immediately when the cavity voltage reaches to the inception voltage due to the availability of free electrons 
for a PD to occur. This process is repeated in the two half cycles of the applied voltage and continues untill 
stopping the applied voltage. Therefore, increasing the cavity diameter increases the PD repetition rates of 
occurrence and increases the maximum PD magnitude. Also, it has been observed that with increasing cavity 
diameter, PD occurs at earlier time and this returns to the availability of electrons from the cavity surface at 
earlier time. 
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Figure 8. Voltage at cavity and PD magnitude with phase angle for one-cycle at 2 mm cavity diameter 
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Figure 9. Voltage at cavity and PD magnitude with phase angle for one-cycle at 4 mm cavity diameter 
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Figure 10. Voltage at cavity and PD magnitude with phase angle for one-cycle at 6 mm cavity 
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Figure 11. Voltage at cavity and PD magnitude with phase angle for one-cycle at 8 mm cavity diameter 


5. COMPARISON BETWEEN MEASURED AND SIMULATED RESULTS 

Table 2 summarizes the measured and simulated results obtained at 30 kV applied voltage, 50 Hz. 
The comparison is based on the maximum PD magnitude in all cavity diameters. It has clearly indicated from 
Table 2 and Figure 12 that the results obtained from simulation model are compatible, to a great extent, with 


those acquired from the experimental measurement for the four studied cavity diameters with a maximum error 
of 8.20 %. 
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Table 2. Measured and simulated results for all cavity diameters 
Cavity diameter (mm) | Maximum PD measured (pC) | Maximum PD simulated (pC) Error + (%) 


2 717.4 658.6 8.20 
4 1295 1236 4.56 
6 3163 3213 1.58 
8 3984 4074 2.26 
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Figure 12. Experimental and simulation results at 30 kV and 50 Hz 


6. CONCLUSION 

In this current work, PD measurement and simulation were performed on a rubber sample with an 
artificial variable cavity diameter to discuss the influence of cavity diameter on the behavior of PD activities. 
PD measurement was carried out by using HFCT technique indicating that the PD magnitude is strongly 
depending on the spatial geometry of the cavity inside the insulating material as well as increasing the cavity 
diameter increases the PD magnitude depending on the applied voltage. Also, a PD simulated model is 
developed using COMSOL Miultiphysics interfaced with MATLAB software for extending a deeper 
understanding of PD occurrence inside a cavity. The simulated model studies the distribution of the electric 
potential and electric field inside the cavity in case of before and after occurrence of PD. Therefore, simulation 
results indicate the PD magnitude and PD repetition rate increase with increasing cavity diameters and PD 
occurs at earlier time for large cavity diameter due to the availability of free electrons emitted from the cavity 
surface at earlier time. The results acquired from the simulation model matches with those obtained from the 
experimental study with a maximum error of 8.20 % based on the maximum PD magnitudes at 30 kV and 50 Hz. 
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